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We study the problem of three ultracold fermions in different hyperfine states loaded into a 
lattice with spatial dimension D = 1,2. We consider S!7(3)-symmetric attractive interactions and 
also eventually include a three-body constraint, which mimics the effect of three-body losses in the 
strong-loss regime. We combine exact diagonalization with the Lanczos algorithm, and evaluate 
both the eigenvalues and the eigenstates of the problem. In D = 1, we find that the ground state 
is always a three-body bound state (trion) for arbitrarily small interaction, while in D — 2, due 
to the stronger influence of finite-size effects, we are not able to provide conclusive evidence of 
the existence of a finite threshold for trion formation. Our data are however compatible with a 
threshold value which vanishes logarithmically with the size of the system. Moreover we are able 
to identify the presence of a fine structure inside the spectrum, which is associated with off-site 
trionic states. The characterization of these states shows that only the long-distance behavior of 
the eigenstate wavefunctions provides clear-cut signatures about the nature of bound states and 
that onsite observables are not enough to discriminate between them. The inclusion of a three-body 
constraint due to losses promotes these off-site trions to the role of lowest energy states, at least in 
the strong-coupling regime. 

PACS numbers: 37.10.Jk, 67.85. Pq, 67.85.-d, 67.80.kb 



I. INTRODUCTION 

The existence of a threshold for bound-state formation 
in the two-body problem is one of the basic differences 
between quantum and classical mechanics^. It is also well 
known that the existence of this finite threshold strongly 
depends on the spatial dimension and on the presence 
(and quantum statistics) of other particles^. 

Even more striking consequences arise however when 
we consider the three-body problem. In this case, ac- 
cording to quantum mechanics, three-body bound states 
are favored in suitable parameter regimes with respect to 
two-body bound states. 

This interplay between two-body and three-body 
physics leads to the so-called Efimov effeclP. Although its 
theoretical prediction by V. Efimov dates back to 1970, 
the incredible experimental advances in ultracold gases, 
which allow to observe and study this physical scenario 
in dilute atomic cloudiP^, have triggered in recent years 
a strong revival of interestP. 

In this paper we address the problem of three-body 
bound state formation in a lattice for three fermions with 
different values of an internal degree of freedom. Our 
calculations show that three fermions can form different 
bound states (see Fig. [lj. This is directly relevant for 
ultracold atomic gases, since the problem under investi- 
gation corresponds to loading three fermions of the same 
atomic species (e.g. 6 Li [THS] or 173 Yb [10]) in different 
hyperfine states into an optical lattice. Similar systems 
have already been studied in current experiments with- 
out superimposing an optical lattice, though the exper- 
imental observation is hampered by three-body losses^. 
Remarkably, by loading the gas into an optical lattice the 



wvww 
wwvw 
vwvwv 
wvww 

FIG. 1: (Color online) Sketch of the particle configurations 
found in the system: (a) onsite trion, (b) off-site trion, (c) 
dimer + atom and (d) unbound atoms. 



effect of the losses could be suppressed, as a large rate 
of on-site three-body loss can prevent coherent tunneling 
processes from populating any site with three particles^] 
For three-body loss, this mechanism wo uld p rovide an ef- 
fective three-body hard-core constraint 12 * 13 !, which sup- 
presses the actual loss events. 

The analogous problem in the absence of a lattice has 
been investigated using the functional renormalization 
group in Ref. [14] , where the authors found a very inter- 
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FIG. 2: (Color online) Ground-state probability distribution uio{xi,X2) for a linear chain of N a = 81 lattice sites for different 
values of the attraction (a) U/t = —0.1, (b) U/t — — 1, and (c) (7/t = —5. 



esting evolution of the system properties for increasing 
attraction. Indeed, in the three-dimensional case, three- 
body bound states appear for smaller values of the at- 
traction than two-body bound states. Therefore there is 
a range of parameters where the ground state is a three- 
body bound state and no two-body bound states are al- 
lowed. In the strong-coupling regime instead a two-body 
bound state turns out to be the stable solution. 

This peculiar phenomenon, which is consistent with 
the Efimov picture, is clearly mirrored in the phase 
diagram of the corresponding many-body proble m 1 14 ^ 15 !, 
though the presence of other fermions triggers pair for- 
mation and superfluidity also at weak-coupling. 

The presence of a lattice is known to induce important 
differences in the many-body phase diagram, especially 
at strong-coupling and in presence of three-body losses. 
Indeed, although the overall picture of the competition 
between superfluid and trimer formation (trionic phase) 
is similar to the continuum case, the latter is stable in 
the lattice at strong-couplingp^lSS unless the three-body 
loss rate 73 is very larg o 11 ! 2 ^. In the continuum instead 
the trionic phas e is b ound to disappear for large values 
of the attractio rPEI. 

In order to better understand this point and the dif- 
ferences to the continuum case, in this paper we take 
a step back and focus on the three-body version of 
the problem investigated in Ref. (2HI- Specifically, 
we consider the case where all the interactions be- 
tween the species are symmetric in a one-dimensional 
chain (ID) and a two-dimensional square lattice (2D). 
We address this problem using a numerical approach, 
combining exact-diagonalization (ED) with the Lanczos 
algorithm 25 , which allows us to obtain information about 
the low-energy sector of the spectrum for much larger 
system sizes than the full ED approach. Moreover we 
also study the effect of including an onsite three-body 
hard-core constraint, which mimicks the effect of a very 
large three-body loss rate 73 in presence of a lattice, as 
explained in Ref. 

The paper is organized as follows. In the next sec- 
tion, we introduce the model and give details about the 



methods. Then, in Section |HI| we discuss the results for 
a one-dimensional chain, first considering the case with- 
out losses, then introducing the three-body constraint 
induced by losses. In Section IV we consider the same 



problem in the two-dimensional case, in order to under- 
stand what is the role of dimensionality in the formation 
of three-body bound states and in the competition be- 
tween two- and three-body bound states. Conclusions 
are drawn in Section Ivl 



II. MODEL AND METHOD 

The system under investigation is described by the fol- 
lowing single-band Hubbard Hamiltonian: 



u 
"2 



(i) 

where (ij) denotes a summation over neighboring lattice 
sites, c i(7 is the fermionic creation operator with hyperfine 

spin a = 1, 2,3 at site i, while fii a = cj CT c io . is the corre- 
sponding number operator, t is the hopping amplitude 
between nearest neighbor sites and assumed as energy 
unit, while U < is the on-site attractive interaction. 

In the context of ultracold gases this Hamiltonian can 
be realized by loading three fermions in different hy- 
perfine states into an optical lattice. Under suitable 
conditions^!, the effect of higher lattice bands can be 
neglected. The case in which the interactions are sym- 
metric corresponds to equal scattering lengths between 
different hyperfine states. In practice this can be realized 
by different hyperfine states of 6 Li in the presence of a 
strong magnetic field, or by alkaline-earth atoms such as 
173 Yb. The Hamiltonian is then (approximately) invari- 
ant under SU (3) rotations in the space of the hyperfine 
states. 

The three-body interaction term with V — > 00 (hard- 
core constraint) is introduced to take the effects of three- 
body losses in the strong loss regime into account ac- 



3 




-10 -8 -6 -4 -2 



u/t 

(a) 

-5 1 ■ ■ .— 




-3 -2.5 -2 -1.5 -1 -0.5 



U/t 

(b) 

FIG. 3: (Color online) Energy spectrum in ID as a function of 
U/t. In Panel (a) we plot the full spectrum obtained within 
ED for N 3 = 19 sites, while in Panel (b) we plot the first 
N = 71 low-lying eigenvalues obtained within the Lanczos 
algorithm for N s — 59. Different bands are plotted with dif- 
ferent color (the trion band is green, the off-site trion band 
is blue, the dimer band is red and the free particle band is 
purple). For weak interactions bands are not separated and 
therefore colors are only formal. 

cording to Refs. |11I13| . On the other hand V = 
corresponds to the case when three-body losses are neg- 
ligible. A finite value of V however, is not equivalent 
to the regime of moderate loss, where actual losses can- 
not be avoided, and the effective Hamiltonian description 
above fails. Therefore we do not consider this situation 
here. 

By using an exact diagonalization approach we were 
able to calculate the full spectrum of the problem for 
system with up to N s — 21 lattice sites for the linear 
chain in ID and A^ s = 16 (4 x 4) lattice sites for a two- 
dimensional square lattice. Since, however, the size of 
the system addressed is relatively small (especially for 
the 2D case where the finite-size effects are expected to 



be more severe), we also used the Lanczos algorithm with 
full reorthogonalization to calculate low-lying eigenvalues 
and eigenvectors for system sizes up to N s — 169 (13 x 13) 
lattice sites. 

We also carefully studied the effect of closed and peri- 
odic boundary conditions. All the results shown in this 
paper are obtained with periodic boundary conditions, 
which preserve the translational invariance of the lattice 
and minimize finite-size effects. 



III. RESULTS: ONE-DIMENSIONAL LATTICE 

In this section we present numerical results for a one- 
dimensional chain. 

For the n-th eigenvalue of the system E n , the corre- 
sponding eigenstate can be expressed as 

IV'n) = ^n(xi,X 2 ,X 3 )cl lA cl^ 2 cl 3}3 \0) (2) 

Xx.X2.X3 

where the same notation as above is used. 
\ , \l) n {xi,X2,x 3 )\ 2 provides us with the probability of 
having the particles with the hyperfine-state a = 1,2,3 
at the positions x a in the n-th energy level. Since the 
system is translationally invariant, we can consider the 
following 3-particle probability distribution: 

oj n (x 1 ,x 2 ) = ^ \^ri(xi + x 3 ,x 2 + x 3 ,x 3 )\ 2 (3) 

which is only a function of the relative positions of the 
fermions with the hyperfine-state 1 and 2 with respect 
to the fermion with hyperfine-state 3. We can anticipate 
that our results will prove that only the long-distance be- 
havior of the many-body wave-function (and therefore of 
uS) provides clear-cut signatures of bound-states forma- 
tion. This is a crucial difference between the results pre- 
sented in the next subsection and Ref. [24], where only 
the energy spectrum and onsite observables were investi- 
gated, and also the main reason we arrive at somewhat 
different conclusions. 



A. Lossless regime, V — 

We start by considering the case where three-body 
losses do not play any role, i.e. we take V = in Eq. 
0. 

First we study the evolution of the ground state prop- 
erties of the system as a function of the interaction. 
As shown in Fig. [2] for all values of the interaction 
the ground-state probability distribution u>o(xi,X2) has 
a pronounced maximum at X\ — x 2 — and a fast de- 
cay with the distance away from this point. We also 
found that, already for a relatively small system size, the 
tails of the probability distribution at large distances are 
essentially unaffected by the size of the lattice. Their 
fast decay is consistent with the formation of three-body 
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FIG. 4: (Color online) Probability distribution u} n {x\,x%) corresponding to several excited eigenstates for N s — 19 and U/t = 
—20: (a) excited on-site trion, (b) lowest energy off-site trion, (c) highest energy off-site trion from the upper off-site trionic 
sub-band, (d) lowest energy and (e) excited dimer states, and (f) unbound particles state. The large value of the attraction is 
chosen to be deep in the strong-coupling regime where the band structures are well separated. 



bound states (trions) for arbitrarily small interactions. 
Therefore, similar to the situation for for two-body bound 
states in continuum space^, there is no threshold for 
three-body bound states formation in ID. It is easy to 
realize that at weak-coupling their wavefunction has a fi- 
nite spread in space over many lattice sites and therefore 
trions can be seen as local objects only for large values 
of the interaction. 



The eigenvalue spectrum shown in Fig. 3(a) allows 
to distinguish three main structures which are clearly 
separated from each other only for strong interactions, 
while they overlap at weak-coupling. We found that 
the lowest band separates from the second band for 
\U\ = U = 1.78t in agreement with Ref. (23]. More- 
over a further study of the low-energy spectrum using the 
Lanczos method (see Fig. |3(b)[ ) shows that this threshold 
converges very quickly as a function of the system size. 

We then characterize the eigenstates through a de- 
tailed study of uj n . As one would expect, the an alysis 
reveals that the states of the lowest band (see Fig. 4(a) ) 
are trionic in nature and oj n is maximal in the center. 
Moreover, since trions are already essentially local ob- 
jects for \U\ = Uq, this provides a simple explanation for 
the observed system size- independence of Uq. 




FIG. 5: (Color online) Probability of triple occupation 
oj n (0,0) as a function of the eigenstate label n for increasing 
energy and different values of the interaction: U — —t blue 
dashed line with squares, U/t = —6 red solid line with circles. 
The data shown correspond to a linear chain of N s = 19 



side the second band 
of the plot in Fig 
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which is not visible on the scale 
Indeed the second band is fur- 



More surprising is the existence of a fine structure in- ther split into three sub-bands which are separated by 
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extremely small gaps at strong-coupling of the order of 
t/2. Despite being so close in energy these sub-bands 
accomodate states which are very different in nature. 

The central sub-band is composed of proper dimer 
states where the three-body bound states are dissociated. 
This is clearly evident in Figs. |4(d)| and |4(e) where u> n 
corresponding to different states in the dimer sub-band 
are shown. In these cases uj n (xi, X2) docs not decay any- 
more at large distances and is maximal along the lines 
xi = 0, x 2 = or x\ = x 2 , being very small in the origin 
xi = X2 = 0. These eigenstates correspond to configu- 
rations where a two-body bound state coexists with an 
unbound particle. 

The states below and above the dimer sub-band show 
instead a completely different behavior, since their prob- 
ability distribution has a fast decay at large distances, 
as shown in Figs. |4(b)| and 4(c) This indicates that 
they are excited three-body bound states. This indicates 
that they are excited three-body bound state, i.e. the 
third particle still stays bound to a dimer, although their 
energies are very close to dimer states. The probability 
distribution at x\ — x-i = is however substantially sup- 
pressed as in the dimer states, while is maximal when a 
pair sits in one of the nearest-neighbors sites of the third 
particle. We therefore refer to these states as off-site tri- 
ons, in order to distinguish them from the states in the 
lowest band. For distinction we will refer from now on to 
the latter as on-site trions. 

Our interpretation is that in the case of off-site trions 
the increase of local quantum fluctuations with increasing 
energy is strong enough to reduce the probability of find- 
ing three fermions on the same site w n (0, 0) basically to 
zero even at strong-coupling, while the three-body bound 
state character, which is related to the large-distance be- 
havior of uj n , is not changed. Indeed the probability dis- 
tribution of dimer states and off-site trions are pretty 




FIG. 6: (Color online) Ground-state probability of triple oc- 
cupancy o;o(0,0) plotted as a function of the interaction for 
different values of the system size: N s = 21,3f,41 (blue 
solid line with squares, red dashed line with circles and green 
dashed-dotted line with triangles respectively.) 
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FIG. 7: (Color online) Dispersion relation for (a) on-site and 
(b) off-site trionic bands for a lattice with N s — 19 and 
U = — 20t. The dispersion relations for both on-site and off- 
site trions are well described by the free particle dispersion 
e(k) = — 2Jeff cos(kl) — fi c ff, where I is the lattice spacing 
assumed as length unit. For on-site trions the results are in 
good agreement with the perturbative result J e g = 3t 3 /2U 2 = 
0.00375* and fj, eS = -SU - 6t 2 /U = 60.15i For off-site 
trions instead J c g oc U~ . 



similar at short distances, and this probably also explains 
why they are so close in energy. We were not able to find 
simple arguments to explain the existence of off-site tri- 
ons and also the size of the gaps between the off-site tri- 
onic and dimer sub-bands, which are for example absent 
in D = 2, as shown in Sec. |IV| 

Finally none of the states of the highest band in Fig. 



3(a) are bound, and the probability distribution w n is es- 
sentially flat and small over the whole system, as evident 
in Fig. |4(f )| We refer to this structure as the unbound 
particles band. 

It is thus interesting to observe that having a maximum 
in the probability of finding three-particles in the same 
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FIG. 8: (Color online) Energy spectrum in ID including the 
three-body constraint (V/t = 10 ) as a function of U/t for 
N 3 — 19. Different bands are plotted with different color (the 
off-site trion band is blue, the dimer band is red and the free 
particle band is purple). For weak interactions bands are not 
separated and therefore colors are only formal. 



lattice site is not essential to three-body bound state for- 
mation, which is instead unambiguously related to the 
long distance behavior of the probability distribution. 

In order to better understand this point we further 
studied w n (0, 0), i.e. the probability to find all three par- 
ticles at the same lattice site. This parameter was used in 
Ref. [23]in order to identify trions and hence called trionic 
weight by the authors. From the considerations above it 
is however clear that the existence of off-site trions and 
the non-local character of the trionic wave-function at 
weak-coupling makes this interpretation inconclusive. 

Indeed, the very irregular behaviour of uj n (0, 0) shown 
in Fig. [5] in the weak-coupling regime is simply due to 
the fact that the onsite trionic band overlaps with the 
higher bands. Therefore some of the states shown in Fig. 
[5]are onsite trions, while others are still trions but off-site 
in character or just dimers. For small U these states are 
very close in energy and only the detailed study of the 
eigenstates that we performed can distinguish between 
them. Even though the states in Fig. [5] with a large 
u! n (0,0) are indeed on-site trions, a decrease in w n (0, 0) 
does not always correspond to an unbound state because 
of the presence of off-site trions. Moreover, in the weak- 
coupling regime there are dimer or even specific unbound 
particle states which have a probability for triple occu- 
pation even larger than for three-body bound states. 

Our point made above is dramatically evident at strong 
coupling, where the onsite trionic band is clearly sepa- 
rated in energy from the higher bands. Indeed, while 
w n (0, 0) is close to one for all the states in the onsite 
trionic band, indicating that onsite trions are essentially 
local objects for this value of the coupling, the same pa- 
rameter drops to a very small value for off-site trions, 
which are still three-body bound states. Within the on- 
site trionic band, we also observe as in Ref. [53] that 
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FIG. 9: (Color online) Energy spectrum in 2D as a function 
of U/t. In panel (a) we plot the full spectrum obtained within 
ED for N s — 16(4 x 4) sites, while in panel (b) we plot the first 
N — 100 low-lying eigenvalues obtained within the Lanczos 
method for N s = 81(9 x 9). Different bands are plotted with 
different color (the trion band is green, the off-site trion band 
is blue, the dimer band is red and the free particle band is 
purple). For weak interactions bands are not separated and 
therefore colors are only formal. 

excited trions have a slightly higher value of w„(0,0), 
compared to the ground state. This is probably related 
to the higher potential energy contribution in the excited 
states and to the different relevance of quantum fluctua- 
tions. 

In conclusions local observables like w„(0, 0) are unable 
to distinguish between an off-site trionic and dimer states 
in all interaction regimes although their wave-functions 
are completely different, as shown above. Therefore 
w„(0, 0) should in general not be taken as an indicator of 
three-body bound state formation. 

As shown in Fig. |6j we found the parameter ujq(0, 0) in 
the ground state to increase smoothly with \U\ indepen- 
dently of the system size in contrast to Ref. [23] , where 
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FIG. 10: (Color online) Ground-state probability distribution 
Qo D (di, da) for a a 2D square lattice with N s = 169(13 x 
13): (a) weak-coupling (U/t = —0.1), and (b) strong-coupling 
(U/t — —4.8) as a function of the Manhattan distances di and 
di (see main text). 



a step-like behavior was observed. 

In the strong-coupling limit we expect that trions be- 
have as composite particles and it makes sense to in- 
vestigate their dispersion relation. The results shown in 
Fig. [7] prove that trionic states are distributed in their 
band as free-particles in the tight-binding approximation. 
For onsitc trions the effective chemical potential /i c g and 
hopping J e fj agree very well with perturbative estimates 
given in Ref. [2D]. For off-site trions instead J e g oc U^ 1 



B. Strong-loss regime, V = oo 

As already mentioned in the introduction, realizations 
of three-component Fermi system in ultracold atomic 
gases are plagued by three-body losses, which strongly 
limit the lifetime for observing the interesting phenomena 



under discussion here. However, as shown in Refs. [22123] 
in the presence of a strong loss rate 73 jt ^> 1 in a lat- 
tice, actual losses are suppressed and the system can be 
described by the Hamiltonian in Eq. [T] provided one 
includes an effective three-body constraint (V — > 00). 
Therefore, in order to understand the effect of a strong 
three-body loss rate on three-body bound state forma- 
tion, we include the additional three-body term in the 
Hamiltonian and set V/t = 10 3 . As proven in Ref. [2D], 
this setup is indeed essentially indistinguishable from the 
real constraint. 

The energy spectrum of the system in the presence of 
the three-body constraint is shown in Fig. [8] Due to the 
constraint, triple occupancies are forbidden and therefore 
onsite trions are suppressed (i.e. their energy is lifted by 
an amount proportional to V). 

One could then wonder if an onsite three-body con- 
straint is sufficient to rule out the possibility of three- 
body bound states in the presence of an attractive two- 
body potential which is also onsite. The previous consid- 
erations on the existence of off-site trions clearly indicate 
that the probability of having triply occupied sites is not 
necessarily correlated with the formation of three-body 
bound states. Our results fully confirm this expectation 
since in the constrained case we found the ground state 
of the system to be always a three-body bound state but 
off-site in nature. As before the off-site trionic states are 
separated by a small energy gap from the dimer states. 
Since both off-site trionic and dimer states have similar 
values of w„(0,0) already in the unconstrained case, the 
energetic hierarchy between them is only weakly affected 
by the constraint, despite of the suppression of triple oc- 
cupancies. The existence of off-site trions and hence of an 
off-site trionic phase has been also observed in Ref. [22], 
where the many-body version of the problem investigated 
here was addressed in one spatial dimension. 




FIG. 11: (Color Online) Average size p of the three parti- 
cle system as a function of U/t for different values of lattice 
size N B = 25, 49, 81, 121, 169. The inset shows the ratio p/L, 
where L is linear size of the system. 
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FIG. 12: (Color online) Probability distribution u)„ (di, d,2) in a 2D square lattice with N s = 81(9 x 9) corresponding to (a) an 
excited on-site trion, (b) an off-site trion, and (c) a dimer. The interaction is again chosen to be deep in the strong-coupling 
regime (U/t = —50) where the different bands are well separated. 



Therefore, at least in ID, three-body bound state for- 
mation is a robust phenomenon which takes place for 
arbitrarily small values of the interaction and is robust 
against three-body losses. Moreover, the presence of off- 
site trions show that quantum fluctuations can strongly 
suppress triple occupancies without necessarily inducing 
the breaking of trions in favor of dimers. In order to un- 
derstand how much these conclusions are related to the 
peculiarity of the one-dimensional case, in the next sec- 
tion we will address the same problem in the context of 
a two-dimensional square lattice. 



IV. RESULTS: TWO-DIMENSIONAL LATTICE 

As in the one-dimensional case, we first investigate the 
case where no losses are present and V = 0. The com- 
plete energy spectrum of the system on a square lattice 
with N s = 16 (4 x 4) lattice site as a function of U is 
shown in Fig. 



9(a 



Despite the same overall struc- 
ture as in the one-dimensional case emerges from the 
picture, finite-size effects arc clearly not-negligible. The 
fast growth of the Hilbert space with the number N s of 
lattice sites strongly constrains the linear size of the lat- 
tice which we can handle within full ED and even within 
the Lanczos algorithm we were only able to address up to 
N s = 169 (13 x 13) sites (see Fig. |9(b)[ ). Since finite-size 
effects are expected to scale with the dimension D as a 
function of N s , it is clear that increasing dimension- 
ality is extremely unfavorable for our numerical approach 
and a careful extrapolation to the thermodynamic limit 
is required for all observables. 

As in ID the lowest band separates from the higher 
bands only beyond a critical coupling. The thresh- 
old value is roughly twice larger than the value in ID 
(\U§ D \/t « 3.4 for N s = 16), suggesting that they could 
converge to the same value when rescaled with the band- 
with W oc D. The prefactor is however slightly bigger 
{\U^ D \/t 3.6) in the data obtained within the Lanczos 



approach on a lattice with N s = 169 shown in Fig. |9(b)[ 
suggesting that there could still be sizable finite-size ef- 
fects affecting the threshold value. 

It is worth pointing out that the lowest band now fol- 
lows rather closely the higher bands in a finite range 
of couplings, as evident in Fig. |9(b)| Since this could 
suggest the existence of a finite threshold for three- 
body bound state formation which is absent in ID, we 
also studied in detail the eigenvectors and the proba- 
bility distribution \tp (x\, 2/1, X2, 2/2, £3, 2/3 ) | 2 in the two- 
dimensional case. 

By using again the translational invariance, we obtain 
the following reduced probability distribution: 

, .2Di 



(an, 2/1, 2:2,2/2) 

^2 \lp n (xi +£3,2/1 



(4) 

+ y 3 ,x 2 + 2:3,2/2 + 2/3, x 3 , 2/3) I 2 ■ 



Since, however, oj„ is still a function of four vari- 
ables, we further reduce them by averaging the proba- 
bility distribution lj^ d (x i, V i, X2, 2/2) over all sites at the 
same Manhattan distance 2 ^ and define 



-2D 



{di,d 2 ) = Vw 2 / 3 ^!, 2/1, 2:2,2/2) 



Ax 1 \ + \y 1 \,d l 0\x 2 \ + \y2\A 2 



CdxCd.2 

(5) 

In the expression above the summation runs over —L/2 < 
2:1,2/1,2:2,2/2 < L/2 and Cd is the number of the lattice 
sites at a Manhattan distance d from the third-species 
fermion. While this allows us to produce meaningful 
three-dimensional plots also in 2D, it clearly averages out 
a part of the information contained in w 2D . 
Moreover, we introduce the following quantity 



P 



= ^ d 1 Cd 1 Cd 2 u>2D{di,d 2 ) 



(6) 



di ,d 2 



which measures the average diameter of the three-particle 
cluster, as clarified below. 

We first consider the evolution of ground state prop- 
erties as a function of the interaction. In the strong- 
coupling regime shown in Fig. |10(b)| where the bands 
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in the spectrum are well separated, the probability dis- 
tribution has a clear maximum at di — di — and a 
very fast decay with distance. In this regime the average 
size of the cluster p plotted in Fig. [TT] is extremely small 
and essentially independent of the system size. Thus we 
can confidently argue that at strong-coupling the ground 
state is an on-site trion. 

Due to the relevance of finite-size effects it is much 
harder to draw analogous conclusions in the weak- 
coupling regime, as shown in I0(a)| Indeed in this case 
the probability distribution stays finite at the edges of the 
lattice even for N s = 169, as it was the case for dimer 
states in ID. On the other hand, since also three-body 
bound states are expected to have a rather large spread 
of the wavefunction in the weak-coupling regime, the fi- 
nite value found at the border is equally likely to be a 
residual finite-size effect. 



c 
LU 




free particle band 



-30 



^dimer band 
\>ff site trion bands 



-20 



-15 



-10 

u/t 



FIG. 13: (Color online) Energy spectrum in 2D including the 
three-body constraint (V/t = 10 3 ) as a function of U/t for 
7V S = 16(4 x 4). Different bands are plotted with different 
color (the off-site trion band is blue, the dimer band is red 
and the free particle band is purple). For weak interactions 
bands are not separated and therefore colors are only formal. 



In order to further clarify this issue we investigated 
the ratio p/L as a function of the interaction strength 
for different values of N s (see inset in Fig. 11). One 



would expect indeed a different scaling with the system 
size for trions and dimer states. 

Our data clearly show that at strong-coupling p/L 
scales to zero as expected for three-body bound states. 
At weak-coupling instead p/L seems at first-sight to con- 
verge to a finite value, with a steep crossover between 
the two regimes around \U\ — U c , where we identify U c 
as the inflection point of p(U) where (d 2 p/dU 2 )N s = 0. 
This threshold value clearly decreases in modulus with 
increasing system size, but we are unable to determine 
whether U c converges to zero or a finite value in the ther- 
modynamic limit, due to the finiteness of our system. 

It is worth mentioning that U C (N S ) seems to be com- 
patible with a logarithmic scaling with the system size, 
as evident from fitting our numerical results. While this 



makes virtually impossible to draw a clear-cut conclusion 
only based on our exact-diagonalization method, it has 
instead a relevance for real systems. Indeed this would 
imply the existence in finite 2D systems, such like e.g. in 
ultracold gases experiments, of a finite threshold for the 
crossover to a three-body bound state for all the practical 
purposes. Similar arguments have already been used in 
the context of ultracold gasea^ to explain the existence 
of quasi-condensates in bosonic systems in 2D, in anal- 
ogy with the existence of a finite magnetization at finite 
temperature in the finite two-dimensional XY modeP^. 

The study of excited states resulted in an outcome very 
similar to the ID case (see Fig. 12 1. The main difference 



compared to the one-dimensional case is that the fine 
structure within the second band in the energy spectrum 
is now absent and off-site trionic and dimer states are 
mixed. This is probably due to the larger bandwidth 
in 2D than in the one-dimensional case, which causes 
the different sub-bands to overlap. The lowest band still 
accomodates only on-site trions (see Fig. 12(a) ) and the 
lowest state in the second band is still an off-site trion (see 
Fig. 12(b) ), at least in the strong-coupling regime where 



the nature of the three-body bound state is evident in 
our data. 

Finally we set V/t — 10 3 and study the effect of the 
loss-induced constraint on the two-dimensional system. 
The picture which emerges (see Fig. 13 ) is pretty similar 



to the one-dimensional case and the lowest-energy off-site 
trion is promoted to the role of ground-state in presence 
of constraint. Due to the greater delocalization of the off- 
site trionic wavefunction, it is even harder to assess the 
formation of three-body bound states at weak coupling 
in this case. 



V. CONCLUSION 

In this work we addressed the problem of three 
fermions in different internal states loaded in a lattice, 
using an exact diagonalization approach. We studied 
the attractive and full S'[/(3)-symmetric case, focusing 
on the competition between three-body and two-body 
bound state formation. 

Our results show that in one dimension three-body 
bound states (trions) are formed for arbitrarily small at- 
traction between the fermions and they remain the low- 
est energy states for every value of the interaction. We 
also observed the existence of excited three-body bound 
states which share the same short-distance structure and 
typical energy of two-body bound states (dimers), while 
they still feature the characteristic long-distance decay in 
the wavefunction of three-body bound states. The prob- 
ability of triple occupancy is strongly suppressed in these 
states, in contrast to the "canonical" trionic states. We 
therefore refer to them as off-site trions in order to dis- 
tinguish them from the states in the lowest energy band 
(on-site trions). 

These results prove unambiguously that the knowledge 
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of the spectrum and onsite observables alone, like in Ref. 
|24j . do not provide a clear-cut distinction among two- 
body and three-body bound-states and can thus lead to 
misleading conclusions. In this respect only the detailed 
study of the long-distance behavior of the 3-body eigen- 
states is able to provide a proper characterization of the 
problem. 

The introduction of a three-body constraint mimicking 
three-body losses in the strong-loss regime as in Ref. p!] ; 
inhibits on-site trions formation and promotes off-site tri- 
ons to the role of lowest energy states. 

In two dimensions, there is evidence in our data of a 
logarithmic scaling with the system size of the thresh- 
old interaction for three-body bound states formation. 
The greater relevance of finite-size effects in the two- 
dimensional case did not allow us to assess if there is a 
finite threshold in the thermodynamic limit, although the 
scaling we found implies the existence of a crossover at 
a finite value of the interaction for all practical purposes 
in finite systems, like e.g. in ultracold gases experiments. 



This also suggests that a finite threshold is likely to exist 
in the three-dimensional case, which we cannot address 
at present using the methods in this paper and which we 
plan to address in a future work by using Monte Carlo 
techniques. 
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